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Abstract 

For fermionic fields on a compact Riemannian manifold with boundary one has a choice 
between local and non-local (spectral) boundary conditions. The one-loop prefactor in the 
Hartle-Hawking amplitude in quantum cosmology can then be studied using the generalized 
Riemann (-function formed from the squared eigenvalues of the four-dimensional fermionic 
operators. For a massless Majorana spin-| field, the spectral conditions involve setting to 
zero half of the fermionic field on the boundary, corresponding to harmonics of the intrinsic 
three-dimensional Dirac operator on the boundary with positive eigenvalues. Remarkably, 
a detailed calculation for the case of a flat background bounded by a three-sphere yields 
the same value £(0) = ^ as was found previously by the authors using local boundary 
conditions. A similar calculation for a spin-f field, working only with physical degrees 
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of freedom (and, hence, excluding gauge and ghost inodes, which contribute to the full 
Becchi-Rouet-Stora-Tyutin-invariant amplitude), again gives a value £(0) = —1|| equal to 
that for the natural local boundary conditions. 


Phys. Rev. D 44, 1713 (1991). 
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I. INTRODUCTION 

In recent work by the authors , 1,2 one-loop effects in quantum cosmology for fermionic 
fields have been studied using local boundary conditions. In the case of a massless spin- -f 

field (vp A , A ^ in a Riemannian background, which we shall refer to loosely as a Majorana 

spin- \ field , 2 the simplest natural local boundary conditions are (using two-component 
spinors) 

e n A i\) A = e ip A , ( 1 . 1 ) 

where e n AA is the spinor version of the unit Euclidean normal e n M to the boundary 
and e = ±1. Note that the primed field ifi A is taken to be independent of i/j a , not re¬ 
lated by any conjugation operation. A first-order differential operator for this Riemannian 
boundary-value problem exists which is symmetric and has self-adjoint extensions. One 
can then study the generalized Riemann zeta-function formed from the squared eigenvalues 
of the Dirac operator. The value of £(0) yields the one-loop divergence of the quantum 
amplitude for the Hartle-Hawking quantum state subject to these boundary conditions. 
Further, £(0) determines the scaling of the one-loop amplitude : in the case of a flat Eu¬ 
clidean four-dimensional background geometry bounded by a three-sphere of radius a, the 
one-loop amplitude scales as a - ^°) for a fermionic field (in the case of a scale-independent 
measure). A direct calculation 2 for a massless Majorana spin- 7 ; field with the boundary 
conditions (1.1) on a three-sphere in flat four-space gave the value £(0) = Such bound¬ 
ary conditions are of interest because they are part of a supersymmetric family of local 
boundary conditions for both fermions and bosons , 1-5 so that one can check whether or not 
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the one-loop divergences in the Hartle-Hawking amplitude cancel in extended supergravity 
theories. 

Because of the first-order nature of the fermionic operators, one has a choice between 
local boundary conditions such as Eq. (1.1) and non-local (spectral) boundary conditions. 
While spectral boundary conditions are not in any obvious way related to supersymmetry, 
they are nevertheless of considerable mathematical interest, and are the subject of this 
paper. Their mathematical foundations lie in the theory of elliptic equations and in the 
index theory for the Dirac operator. 6 To illustrate these boundary conditions, consider 

again the case of a massless Majorana spin-i field (^ip A , A j in the region of flat Euclidean 
four-space bounded by a three-sphere of radius a. Denote by r the Euclidean distance from 
the centre of the sphere. Then the field (^> A , L A j may be expanded in terms of harmonics 
on the family of spheres centred on the origin, 2,7 as 


. 3 

A = 4 — a ™ m np(r)p nqA + r np {r)a nqA 


2tt 


( 1 - 2 ) 


npq 


^ 2tt ^ n 

npq 


+ r np ( T )a nqA ' 


(1,3) 


In the summations, n runs from 0 to oo, p and q from 1 to (n + 1)(n + 2). The a^ q are a 
collection of matrices introduced for convenience, where, for each n, cd’' 1 is block-diagonal 


in the indices pq , with blocks y | J harmonics p nqA have positive eigenvalues 

\(n+ |) of the intrinsic three-dimensional Dirac operator e nAA'G BA Dj on the three- 
sphere, while the harmonics a nqA have negative eigenvalues — |(n + |). Here e BA ^ is 
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the spinor version of the orthonormal spatial triad on the three-sphere, and ^Dj is the 
three-dimensional covariant derivative (j = 1,2,3).' Similarly, the harmonics a nqA have 
positive eigenvalues \{n -\- |) of the corresponding three-dimensional operator on primed 
spinors, and the harmonics ~p nqA have negative eigenvalues —^(n+ |). This expansion 
can be summarized more simply as 


= 4>f+) + yl) 


(1.4) 



yl) + y+, 


(1.5) 


where the (±) parts correspond to positive and negative eigenvalues, respectively, for the 
intrinsic three-dimensional Dirac operator. 

In studying the classical boundary-value problem for the massless Dirac equation, one 
finds that classical solutions corresponding to boundary data with a non-zero coefficient 
r np (a) or m np (a ) diverge as a negative power of r near the origin.' Boundary data with 
a non-zero coefficient of m np (a ) or r np (a ) yield a regular solution of the massless Dirac 
equation, proportional to a positive power of r. Thus the classical boundary-value problem 
is only well-posed if one specifies the m np (a ) and r np (a), but not the remaining data. In the 
case of a general manifold with boundary, knowledge of the spectrum of the intrinsic three- 
dimensional Dirac operator is necessary if one wishes to compute the ^-invariant which 
gives a boundary contribution to the index of the Dirac operator for the manifold with 
boundary. 6 [In the generic case, the index will be non-zero, and the classical boundary- 
value problem will not be well-posed. 6 ] In this paper we are however not concerned with 
the signature, but rather, as explained above, with the zeta-function £(s) formed from the 
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eigenvalues of the four-dimensional Dirac operator, subject now to boundary conditions in 
which and ^ A + ) are specified on the boundary in our flat-space example. 

Thus, just as one has a well-posed classical problem with these boundary data, one 
similarly expects that the analogous quantum amplitude, the Hartle-Hawking path integral 

K HH = J e~ lE D^ A D^ a ' (1.6) 

for the fermions, is naturally studied by taking spectral boundary conditions in which 
and ip A + x, or equivalently the m np {a ) and r np {a) in our example, are specified on the 
boundary. Here the Euclidean action is 

Ie = \ f d^xy/g [V {Vaa^ A ) ~ (Vaa'V) V\ + I B • ( 1 - 7 ) 

The fermionic fields are taken to be anti-commuting, and Berezin integration is being used. 7 
With our conventions, the Infeld-van der Waerden connection symbols a AA are taken to 
be cto = — " 75 ^’ ai ~ 'tI ~ 1,2,3), where are the Pauli matrices. A boundary term 
7 b, discussed in Ref. 7, is needed in general. In our simple example of the three-sphere, 
the Hartle-Hawking amplitude is then a function Khh Wn P (a), Di P (o.) of the spectral 
boundary data. The one-loop properties of this amplitude can be studied without loss of 
generality by setting the allowed boundary data m np (a), r np (a ) to zero, Vn,p, so that the 
boundary conditions become 

<+)=0 . 7 ( 7 , =0 ■ ( 1 . 8 ) 
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[A similar simplification was made in Ref. 2 in choosing the local boundary conditions 
(1.1), rather than specifying a non-zero spinor field \/2 e n/ r ip A — e ip A on the boundary.] 
The boundary term //> is zero in this case. 

In Sec. II the action (1.7) is expanded in harmonics, subject to the spectral boundary 
conditions (1.8) on the three-sphere of radius a bounding a region of flat four-space. The 
eigenvalue equation arising in the evaluation of the one-loop functional determinant is 
derived. The resulting zeta-function formed from the squared eigenvalues is related to the 
heat kernel G(T ) for the Laplacian operator on spinors. In Sec. Ill the Laplace transforms 
of the corresponding Green’s functions are derived. The detailed calculations leading to 
the asymptotic expansion of G(T ) as T —> 0 + are described in Sec. IV. These lead to the 
value C(0) = which remarkably is the same as that found previously 2 for a massless 
Majorana spin-1 field using the local boundary conditions (1.1) on the same manifold with 
boundary. An analogous calculation for the spin-| field is sketched in Sec. V, working 
only with the physical degrees of freedom in a particular gauge. This of course excludes 
the contribution of gauge and ghost modes which should appear in the full Becchi-Rouet- 
Stora-Tyutin-invariant path integral. Nevertheless, it is striking that the value £(0) = —1|| 
obtained is again identical to that found using the natural local boundary conditions for 
the physical degrees of freedom. 1 Some comments are included in Sec. VI. 
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II. EIGENVALUES FOR SPECTRAL BOUNDARY 
CONDITIONS ON S 3 

In the case of a region of four-dimensional flat Euclidean space bounded by a three- 
sphere of radius a, we decompose the massless spin-i field (^p A , '4 >A ^ as in Eqs. (1.2,3) 
and impose the spectral boundary conditions (1.8), so that 

m np (a) — 0 , r np (a) — 0 , Vra,p . (2.1) 

The Hartle-Hawking path integral (1.6), with the Euclidean action Ie given by Eq. (1.7), 
can then be studied, with m np (r), r np (r) constrained by Eq. (2.1) at the boundary, but 

m np (a) and r np (a) unconstrained. The physical fields (^> A , V ’ A ) summed over in the 

3 

path integral should at least be bounded near the origin r = 0. Because of the factor r~ 2 
in Eqs. (1.2,3), this implies that 

m np { 0) = r np (0) = m np (0) = r np (0) = 0 , Wn,p . (2.2) 

The action Ie can then be expanded out in terms of harmonics, by analogy with the 
treatment of Ref. 7, as 

00 (»i+l)(n+2) 

Un(m np , m np ) “t - 7 n (r np , r np ) , (2-3) 

n=0 p= 1 

where 

(2.4) 
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and an overdot denotes jjr. Note that, as remarked in the Introduction, no boundary 
term /g of the type described in Ref. 7 appears in the action, because of the boundary 
conditions (2.1,2). 

Because of the degeneracy (n + l)(n + 2) in the label p. and because of the splitting 
of the action Ie in Eq. (2.3) into identical pieces involving ( m np ,fh np ) and (r np ,r np ), the 
complete path integral (1.6) splits into a product of Berezin integrals : 

00 ( r . 2(n+l)(n+2) 

Khh — n{ / d[x]d[x\exp[-I n (x,x)]j . (2.5) 

n=0 ^ 

The boundary conditions in each integration are then 

z(0) = x(0) = 0 , x(a) = 0 , (2.6) 


following Eqs. (2.1,2). 

Equivalently, one can follow the procedure of Ref. 2 and study the eigenvalue equations 

V AA'^Prn j ^AA'^rn ^m^PmA > (2-7) 

naturally arising from variation of the action (1.7), subject to the boundary conditions 
(1.8). The eigenfunctions V’m ) are clearly found by separation of variables, being 

of the form 

t A = rix nfc (r)^ , =r~ix nk (r)p npA ' , (2.8) 

or 

V’n = T~?x nk {T)a npA , i> A ' = T~^x nk {T)a npA ' . (2.9) 
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Here the pair yx n k{r), x n k{r)J is an eigenvector corresponding to variation of the action 
I n (x,x ) of Eqs. (2.4,5), obeying 



Xnk LJ n kX n k 




( 2 . 10 , 11 ) 


subject to the boundary conditions x n k{0 ) = x n k{0 ) = 0, x n k{a ) = 0, where v = 
and E n k = iX n k■ For each n. the index k labels the countable set of corresponding eigen¬ 
functions and -values. 

The Gaussian fermionic path integral (2.5) is then formally proportional to the product 

1 -r / A \ 2(n+l)(n+2) _ 

of the eigenvalues [[ n k J . Here the constant // with dimensions of mass 

has been introduced in order to make the product dimensionless. 2 As shown below, the 
X n j : are purely imaginary, or equivalently the E n k are real. Further, the E n k occur in 

equal and opposite pairs, since if the pair ^x n fc(r), x n k{r )j corresponds to an eigenvalue 
E n k in Eq. (2.10), then clearly the pair (^x n k{x), — x n k(r )j corresponds to the eigenvalue 

/ | \ 2(n+l)(n+2) 

— E n k■ Hence the path integral can also be written as k ( ~ j . The Berezin 

integration rules imply that one should only include those values of k, say k — 1 , 2 ,..., which 
correspond to positive values of E n k. This formal expression must then be regularized using 
zeta-function methods. 

The coupled first-order equations (2.10,11) lead to the second-order equation 



[(n+ 1) 2 ~|] 

r2 



Xnk — 0 


( 2 . 12 ) 
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together with a corresponding equation for x n k{r). The solutions (x n j t(r), x n k{r obey¬ 
ing the boundary conditions (2.6) are 

Xnk — A n k\/x Jn+l(Enk'r') : ( 2 . 13 ) 

%nk A n k\/x Jn+2(dS n k'T~') , (2.14) 

subject to the eigenvalue condition 

Jn+i{E n ko) — 0 , n = 0,1,2,.... (2.15) 

The quantities E^ k are clearly real and positive since they are eigenvalues for the self- 
adjoint problem (2.12), subject to boundary conditions x n k( 0) = x n k{a ) = 0, and since 
Eq. (2.12) involves a positive operator. For a given n, the coefficients A n k in Eqs. (2.13,14) 
can be chosen such that the eigenfunctions x n k(r) are orthonormal in the inner product 
(u, v) = Jq dr as are the corresponding x n k{r). For each n. the action I n (x, x) 

of Eq. (2.4) becomes a diagonal sum over the eigenfunctions x n k , T n fc for k = 1,2,.... 
Performing the Berezin integrations in Eq. (2.5), one arrives at the formal expression 
quoted in the previous paragraph : 

oo oo , „ \ 2(n+l)(n+2) 

=n n m . (*■>«) 

n=0k=l \ b / 

where the E n k (k = 1,2,...) are the positive eigenvalues obeying Eq. (2.15). 
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Following the standard procedure, as, for example, in Ref. 2, the formally diver¬ 
gent infinite product (2.16) is regularized by studying the zeta-function for the squared 
eigenvalues 

OO OO 

■ (2-17) 

n=0 k=1 

Here the degeneracy c4(n) = 2(n + l)(n + 2) is in fact independent of k. The series (2.17) 
converges for Re(s) > 2, and can be analytically continued to a meromorphic function 
with poles only at s — |,1,|,2. The formal expression log (Khh), with Khh given by 
Eq. (2.16), is then evaluated as — |e'(0) — C(0) log/r. 

III. GREEN’S FUNCTIONS AND THE HEAT KERNEL 


The quantity C(0), which gives the divergence and scaling properties of the one-loop 
amplitude, is evaluated by studying the heat kernel, defined for T > 0 by 

OO OO 

cm = E E 2{n+ l)(n + 2)e“ (E " fc)2T . (3.1) 

n=0 k =1 


This is related to the zeta-function £(s) of Eq. (2.17) by 


COO = 


r(E 


poo 

/ dT T S ~ 1 G(T) 

Jo 


(3.2) 


for Re(s) > 2. The heat kernel G(T) has the standard asymptotic expansion 


G(T ) ~ V BiTi 


(3.3) 


i =0 
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as T —*► 0 + , where in particular I3\ — C(0). 8 

In the present case of spectral boundary conditions, the eigenvalue condition 


Jn+i (L/ n ktt') 0 , n 0,1,2,... 


[Eq. (2.15)], with degeneracy 2(n + 1 )(n + 2), is considerably simpler than the eigenvalue 
condition [J n+ i(E n ka)] 2 — [J n+ 2 {E n ka)] 2 = 0, n = 0,1,2,..., with degeneracy (n+l)(n+2), 
found in Ref. 2 for a spin-| field with local boundary conditions on S 3 . As a consequence, 
we can use a more straightforward treatment, following (among others) Schleich 8 and 
Stewartson and Waechter, 9 which involves the Green’s function for the heat equation for 
each n. The value £(0) = ^ which results is, perhaps surprisingly, equal to that found in 
Ref. 2 for local boundary conditions. 

One proceeds by considering, for each n — 0,1,2,..., the Green’s function defined for 
T > 0 by 

OO 

G n (r, t', T) = y ^x nk (T)x nk (r')e~^ Enk ^ T , (3.4) 

fc=l 

with G n (r, r', T) — 0 for T < 0. Here the x n k(r ) are the eigenfunctions of Eq. (2.12), obey¬ 
ing Eq. (2.12) and x n k(0 ) = x n k{a) = 0, and normalized according to f ( Q dr x n k(T)x n i(T) = 
dki-, as described in Sec. II. Here G n (r, r',T) is the Green’s function for the heat equation 


d 

df 


d 2 ((n+l) 2 -|) ~ 

dr 2 t 2 


G n (r , r', T) = 6(r 


t')S(T) 


(3.5) 


It obeys the boundary conditions 


G n (a, r', T) 


G n {r , a, T) = G n (0,r',T) 


G n {ji 0, T) = 0 


(3.6) 
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By setting t — r' and integrating, one recovers the contribution 


Gn(T)= [ dr G n (T,T,T) = f]e 
Jo k =i 


~{E nk ) 2 T 


to the heat kernel 


G(T) = £ 2(n + 1)( 77 . H— 2)G n (T) 


The Laplace transform of the Green’s function, 


G n (r, t', a 2 ) = / dT e a r G n (r, r', T) , 


obeys the differential equation 


d 2 2 ((n + l) 2 — \) 


— a — 


G n (T,r'.a 2 ) = -S(t - t) 


(3.10) 


Following Eq. (3.6), G n (r , r', a 2 ) is zero whenever either r or r' is 0 or a. It can be found 
explicitly in terms of modified Bessel functions (cf Refs. 8-11), as 
O 1 _ > _2\ /_ \h^n+ i( crr <) 


Gn(r, r', a ) = (r<) 2 (r>)2 


4i+i (era) 


/ n+ i(aa)Ji n+ i(ar>) -/ n+ i(crr>)/i n+ i(cra) , (3.11) 


where r> (r<) is the larger (smaller) of r and rh This gives a splitting of the Laplace 
transform G{t,t ', cr 2 ) of the function G(r, r',T), defined as 


G(t, t', T) = £ 2(n + l)(n + 2)G„(r, r, T) , 


(3.12) 
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in the form 

G(t, t', a 2 ) — G f (t, t\ a 2 ) + G 1 (r, t', a 2 ) . (3.13) 

Here 

OO 

G F (t, t 1 , a 2 ) = ^ j 2(n+l)(n + 2)(T < )^(T > )^I n+1 {<7T < )K n+1 (<jT > ) (3.14) 

n=0 

is the “free contribution”, which corresponds to the boundary conditions of vanishing at 
the origin and at infinity. The “interacting contribution” is 

OO jy- / \ 

G 1 (t, t', a 2 ) = -^2(n+ l)(n + 2 )(r<) 5 (r>) 5 ^htL^/ n+1 ( ar< )/ n+1 ( (Tr> ) . (3.15) 

I n +i(aa) 

By studying the larg e-a 2 behaviour of these functions, or the corresponding small -T 
behaviour of G(T) as in Eq. (3.3), one finds £(0). 


IV. DETAILED CALCULATION OF THE INFINITE SUMS FOR SPIN \ 


Following the calculation of Schleich, 8 taking the inverse Laplace transform term by 
term and integrating with respect to r as in Eqs. (3.7,8), the free part of the heat kernel 
is found to be 


G f {T) 



OO 

J2[n(n + 1) I n (y)t 


n= 1 


dy 


(4.1) 
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Using the integral representation of the Bessel functions, 11 one obtains the identity 


OO 

^,nU n (y) = \e y 


Z- 

n= 1 


(4.2) 


which implies 


^ r%r tt 4 

n 2 I n {y)e y dy = 


tJo 


16 T 2 


(4.3) 


Moreover, using the relations 


li 


Io(y)=h(y ) , nJ„(y) = |(/n-i(y)-/ n+1 (?/)) (4.4,5) 


one has 


a 2 

, ^ OO 

2T ^^ 


o 2 

. OO 

2T - ^ 


n=l 


nl n (y)e y dy= (jn-M ~ J»+i(y)) | e V d y 


n= 1 


2 2 
a a 

1 /* 2T 1 /* 2T 

2 L ye ~‘ Uy) dy + 2 L ye ~"l^ dv 


la 2 ^a 2 . fa 2 

-e 2T /n — 

2 2T V 2T 


a 2 

+ \J o e ~ y ( 2 y - l ) I o{y)dy 


(4.6) 


The problem of computing G F (T) is thus reduced to that of computing the right-hand side 
ofEq. (4.6). Using again the identity/ q(?/) = Ii(y), together with /{(y) = Io(y)— y _1 /i(y), 
one has the relations among indefinite integrals : 


ye %(y) dy = 


ye 


'(Jo + Ji) - yV^io + h) + y 2 e~ y (h + J 0 - ^ 


dy 


= y 2 e 2/ (/ 0 + /i)- / ye y (/ 0 + /i)c/y , 


(4.7) 
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giving 


3 ye y I 0 {y)dy = y 2 e y (I 0 + h)-ye y I 0 + / Joe y dy . (4.8) 


In addition, 


J 0 e y dy = 


ye~ y (h+Io-^+ilo + h) (e~ y - ye~ y ) 


dy 


= ye- y (I 0 + h ) 


(4.9) 


Hence the integral on the right-hand side of Eq. (4.6) is found to be 


e y (2y - 1 )I 0 (y) dy = e y { | y 2 I 0 + ( y 2 + y)h \ - y(l o + h) \ ■ (4.10) 


The relations (4.1)-(4.6) and (4.10) imply that, as T —■> 0 + 


p a 9 a 3 a i 

G f (T) -T“ 2 H- —T~i -= 

16 bvTF 8^y^r 


+ 0(VT) 


(4.11) 


In deriving Eq. (4.11), we have used the following asymptotic relations 11 valid as z —» oo : 


_*. 2 r / \ 2:2 1 .. ( 1 

e z 2 Iq(z ) ~ / _ H-== + O —= 

1 7 \ \fz 


(4.12) 


2rM 3 Jz ( 1 

e z 2 Ii(z) ~ - -== + O —-= 

\/2tt 8 y/2n \ \fz 


(4.13) 


The Laplace transform of the kernel of the interacting part is given by 


G 1 (a 2 ) — —a 2 ^ n 2 f(n; aa) — a 2 ^ n/(n; era) , 


(4.14) 


n= 1 


n= 1 
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where 


f(n ; era) = 


I n (aa)K n (aa ) - I' n (cra)K' n (cra) 


aal n (aa ) 


(4.15) 


is the function defined in Eq. (5.9) of Ref. 8, and in Ref. 9. In fact (see Kennedy 12 ) the 
sums in Eq. (4.14) diverge because of the factors of n 2 and n. This occurs because 

we are attempting to take the Laplace transform of a function G(T) which is singular as 
T —> 0 [see Eq. (3.3)]. This difficulty can be avoided by first computing the sums Y2n=i 
for large a 2 , using the asymptotic expansion of f(n;aa) valid uniformly 8-10 with respect 
to n at large era, and then taking the inverse Laplace transform before taking the limit 
N —> ex). The first series in Eq. (4.14) has already been studied in the case of scalar 
fields. 9 In the case of the second series, the Watson transform used in Refs. 8-9 is a source 
of complications, because nf(n; era) is not an even function of n. Instead, we take the 
inverse Laplace transform of the large-cx 2 expansion of nf(n: era), and compute the sum 
(as an asymptotic series valid as T —> 0 + ) with the help of the Euler-Maclaurin formula. 13 

Setting r = -y === , one has the asymptotic series 8-9 



r (1 — r 2 ) 
2 n 


r 4 (1 — r 2 ) 

~ 2 ~ 


r 3 (1 — r 2 )(l — 12r 2 + 15r 4 ) 

8 n 3 

t r 4 (1 - r 2 )(2 - 53r 2 + 168r 4 - 125r 6 ) 
T ,, . 


(4.16) 


valid as a -MX), uniformly in n. This is derived from the uniform asymptotic expansions 
of I u (z),K u (z),r v {z), and K' v (z ) described in Refs. 10,11. Denoting by L/ the inverse 
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Laplace transform, one thus has the asymptotic series 

n 2 T n 2 T 


Li 


nf(n ; era) 


1 

~ — 7 T ne or 
2a 2 


3 

2T2 „-- 

-n 6 e a 


3^Tra* 


n 2 T 


T -2 

+ s^ ne a 


3 T 2 
4a 6 


n 2 T 

n 3 e a 2 


5 T 3 


n 2 T 


5 — 1 T§ 

+ 16^ e " + 6^^ ne a 


n 2 T 


53 T 2 168 T2 

rrre a + tttz~ a 


n 2 T 


30 v/yra 7 


105 y/rta 9 


2 «2 T 
50 T2 7 -^4- 

-n'e a +... 


189 v^Fa 


li 


(4.17) 


Note that, when each term on the r.h.s. of Eq.(4.17) is summed from n = 1 to N, the 
resulting function of T does not always converge uniformly to the sum ) in a 

neighbourhood T G (0,5) [see Eqs.(4.19)-(4.21)]. Nevertheless, a study of the error terms 
shows that it is valid to take the limit N — > oo as in Eqs.(4.19)-(4.21), and then examine 
the small- T behaviour of the resulting contributions to G(T). 

n 2 T 

In order to compute sums of the type Y4i=i m + Oe a 2 where m = 0,1,2,..., we 

can use the Euler-Maclaurin formula 13 

1 f°° 1 ~ 

-E(0) + F(l) + F( 2) + ... - j F(y)dy = -~B 2 F'{ 0) 

- §r"'(0) - §F""'(0)... , (4.18) 
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for the function F(y) — ye a 2 . In Eq. (4.18), the B t denote the Bernoulli numbers. 
Thus we get 


OO 



n= 1 



a 2 1 T T 2 
2T 12 ~ 120a 2 ~~ 504a 4 + 


(4.19) 


The other sums arising from Eq. (4.17) are obtained by differentiating Eq. (4.19) with 
respect to T. This yields 


£ 


_n 2 T 
n 3 e a 2 


n— 1 


a 


2 T 2 


+ 


1 

120 + 


T 

252 a 2 


+ ... 


(4.20) 


E 


_ _n 2 T 
n 5 e a 2 


n =1 


a 6 1 
T3 “ 252 + 


(4.21) 


and so 011 . 

Combining the resulting contribution to the asymptotic expansion of G 1 (T) as T —> 0 + 
with the other piece, arising from the first term in Eq. (4.14), given in Ref. 9, as well as 
the expansion (4.11) of G F (T), one Ends the asymptotic expansion of the heat kernel : 


G{T ) ~ 


—T~ 2 + a 3 
16 



3 

T 2 + a 




T 


1 

2 + 


11 

360 


+ o(VT ), 


(4.22) 


valid as T —> 0 + . In particular, this yields 


C(0) 


11 

360 


(4.23) 


for the spin-i field with spectral boundary conditions on the sphere. 
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V. CALCULATION OF £(0) FOR THE SPIN-§ FIELD 
WITH SPECTRAL BOUNDARY CONDITIONS 


In this section we sketch the corresponding calculation of £(0) for a linearized spin-| 
field subject to spectral boundary conditions on a 3-sphere of radius a. As in the original 
quantum-gravity calculation of Schleich, 8 we work only with physical degrees of freedom 
by imposing a gauge condition and constraints. This will exclude the contribution of gauge 
and ghost modes which appear in the full BRST-invariant path integral. [For a BRST- 
invariant approach to computing (’(0) for a spin-f field with local boundary conditions, 
see Ref. 14], The value £( 0 ) = — 1 || found here for spectral boundary conditions is again 
identical to that found using the natural local boundary conditions, 1 working with the 
same physical degrees of freedom inside the 3-sphere. 

The spin- 1 field, as appearing, e.g., in N — 1 supergravity, is described by a potential 

1 P A I / S j i 11 the Euclidean regime (/ / = 0,1,2,3). In a Hamiltonian treatment, 15 the 

quantities j are the dynamical variables (i = 1,2,3), while (ip A 0 , tp A 0 j appear 

as Lagrange multipliers. The gravitational field is correspondingly described by the tetrad 
e AA . Taking the geometry to be flat, and x° = r to be the radial distance from the origin, 
while x l (i = 1, 2, 3) are coordinates on the 3-sphere, we impose the gauge conditions 

e AA'' l l jA j = 0 > e AA J ’^ A j=° ■ (S- 1 ) 

If in addition we require that the dynamical variables ^ jj A t , - obey the linearized 

supersymmetry constraint equations 15 on the family of 3-spheres centred on the origin, 
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then the expansion of ^ A i^j in harmonics, 16 analogous to Eqs.(1.2,3) for the spin-^ 

field, takes the simplified form 


.,.A 


(5.2) 



r 


2t r 


npq 


a 


pq 


m 


np 


(t)P 


■nqBA'B' 


+ r np (r)p 


nqBA'B' 


e-BB' 


(5.3) 


In the summations, n runs from 0 to oo, and p and q now run from 1 to (n + l)(n + 4). 
For each n, a IA is a matrix again block-diagonal in the indices pq, with blocks 
Here 

pnqABB' = _ p nq(ABC) ^B' 

where n AA — i e n AA is the Lorentzian normal, 2,7 and the harmonic p n q( ABC ) ^ sym¬ 
metric in its three unprimed indices, has positive eigenvalue \(n+ |) of the appropriate 
3-dimensional Dirac operator. Similarly 



^nqBA'B' = _ a nq(A'B'C') n B ci 


(5.5) 


where the harmonic a nq ^ A B c ), symmetric in its three primed indices, also has positive 
eigenvalue |(n + |). The harmonics jp^ ABB and (3™ qBA B are given similarly in terms 
of harmonics a nq ( ABC '> and ~p n, L A B c ) ; which have negative eigenvalues —|(n + |) of the 
3-dimensional operator. 
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Just as in the spin-| case, the simplest natural spectral boundary conditions for spin 
| are 

m np (a) — 0 , r np (a) = 0 , Vra,p . (5.6) 

The requirement that the physical fields be bounded near the origin r = 0, in the Hart.le- 
Hawking path integral for the linearized spin-| theory, implies that rn np ,r np ,m npi r np 
vanish Vn, p at r = 0. The Euclidean action Ie for the linearized spin-| field, working 
only with the physical degrees of freedom given in Eqs.(5.2,3), is analogous to the spin-| 
expression in Eqs.(2.3,4) : 


oo (n+1) (n+4) 

fe = E E 

n —0 p= 1 


In {j^np 5 TY^np) H"~ ixnp 9 Xnp ) 


(5.7) 


where 


I n (x, x) — dr 

Jo 


1 (n + 2 ) ~ 

- [ xx + xx ) - —xx 

2 V J T 


(5.8) 


The resulting Berezin integration for the Hartle-Hawking amplitude Khh gives the formal 
expression 


k hh =n n 


E 


2(n+l)(n+4) 


nk 


n= 0 k=l 




(5.9) 


The E n k (k — 1, 2,...) are the positive eigenvalues of the system 


d_ 

dr 


& ) %nk — i^nkXnk ? 


dr 


^ ) X n k — EnkXnk 5 


(5.10,11) 
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subject to x nk {0 ) = x nk (0) = x nk (a ) 


0, where z> = 


(n+|) 


The solutions are of the form 



(5.12) 


(5.13) 


with the eigenvalue condition 


Jn+2 



0 


n = 0,1, 2,... 


(5.14) 


and degeneracy 2 (n + l)(n + 4). 

The formal expression (5.9) for Khh is then evaluated by studying the corresponding 
zeta-function and heat kernel, by analogy with Secs.II-IV. The formulae in Sec.Ill can be 
straightforwardly modified to allow for the different degeneracy 2(n+l)(n+4), different dif¬ 
ferential equations (5.10,11), and eigenvalue condition (5.14). In particular, Eqs.(3.14,15) 
should be modified by replacing the factor of (n + 2) by (n + 4), and by changing the order 
n + 1 of each modified Bessel function I n+i or K n+1 to the order n + 2. 

For the sake of brevity, we consider only the constant part B\ = £(0) in the expansion 
(3.3) of the spin-1 heat kernel G(T) as T —■> 0 + . The results of Eqs.(4.2)-(4.13) show that 
the free part G F (T ) of the heat kernel for spin |, given by Eq.(4.1) with the factor (n + 1) 
replaced by (n + 3), gives no contribution to £(0). The only contribution to ^(0) for spin 
| arises from the interacting part 

OO OO 

G 1 (a 2 ) = —a 2 y~^(n 2 — 2)/(n; aa) — a 2 n/(n; aa) , (5.15) 

n=2 n =2 
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where /(n; era) is given by Eq.(4.15). As in Sec.IV, these sums diverge, and the calculation 
should strictly proceed by first computing the sums Yln=i f° r large a 2 , taking the inverse 
Laplace transform, and finally taking the limit N —> oo. The contribution from the first 
term in Eq.(5.15) is found, following Refs.8,9,12, by using a Watson transform : 


-a 2 n 2 - 2 )/(n; a a) 

n =2 



'C'-Q 


dv (is 2 


2 )f(is; era) cot 7rz/ 


(5.16) 


Here the contour of integration C' — Q (compare with Fig.l of Ref.8) encloses all poles 
along the real axis, except for those at v — 0, ±1. Using the uniform asymptotic expansion 
of f(is]cra), and computing the contribution from the poles at is — 0, ±1, one finds from 
the large-cr behaviour of Eq.(5.16) that the first term in Eq.(5.15) contributes — 4|y- to 
C(0). Using the results of Sec.IV, the second term in Eq.(5.15) is found to contribute 

to c(o). 

Combining these results, we find 

C(0) = -H (5.17) 

for the linearized spin-| field subject to spectral boundary conditions on the 3-sphere, 
working only with physical degrees of freedom in the gauge (5.1). Just as in the spin-i 
case of Sec.IV, the value of C(0) is equal to that found previously for the spin-| field subject 
to the natural local boundary conditions. 1,2 Specifically, the value (’(0) = —1|| was also 
found for spin |, working only with physical degrees of freedom subject again to the gauge 
conditions (5.1), with the local boundary conditions 

V2 eWi 4 V A i = e A 'i , (5.18) 
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where e = ±1. These local boundary conditions are part of a locally supersymmetric 
family of boundary conditions for different spins, 1,2,5,14 which includes the local boundary 
conditions (1.1) for spin Of course, as remarked earlier, the results of £(0) calculations 
for gauge fields such as spin | will be modified by the contribution of gauge and ghost 
modes in a complete BRST-invariant calculation. 


VI. COMMENTS 

It is striking that the same value £(0) = is obtained for a massless Majorana 
spin-1 held on a ball in Euclidean 4-space, bounded by a sphere of radius a, whether local 
boundary conditions (1.1) or spectral boundary conditions are imposed. The different 
eigenvalue conditions, respectively, 

[J n+ i(Ea)] 2 — [J n+ 2 (Ea)] 2 — 0 , n = 0,1, 2,..., degeneracy (n + l)(n + 2) , (6.1) 

and 

J n+ i(Ea) — 0 , n = 0,1, 2, degeneracy 2{n + l)(n + 2) , (6.2) 

offer no obvious explanation of this equality- although they suggest an alternative approach 
through studying the asymptotic distribution of eigenvalues. The same holds for the value 
C(0) = —1|| found for the spin-| held (taking only physical degrees of freedom), both for 
local and spectral boundary conditions. There the eigenvalue conditions are, respectively, 

[J n+2 {Ea)] 2 - [J n+3 (Ea)] 2 = 0 , n = 0,1,2,..., degeneracy {n + l)(n + 4) , (6.3) 
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and 

J n+ 2 (Ea) = 0 , n = 0,1, 2, degeneracy 2{n + l)(n + 4) . (6.4) 

The eigenvalue conditions for fermions, subject to spectral boundary conditions on 
the 3-sphere, are in fact similar, but not identical, to the eigenvalue conditions for bosons, 
subject to local (Dirichlet) conditions on the 3-sphere. For example, in the case of a scalar 
field, 9,12 Dirichlet conditions give 

J n+ i(Ea) — 0 , n = 0,1, 2, ..., degeneracy (n + l) 2 . (6.5) 

For a Maxwell field, taking only physical degrees of freedom, 1 ' Dirichlet (magnetic) bound¬ 
ary conditions give 

J n+2 (Ea) = 0 , n = 0,1, 2,..., degeneracy 2{n + l)(n + 3) . (6-6) 

Analogous equations hold in the spin-2 case. 8 Because of the equality of the £(0) values 
for fermions with both local and spectral boundary conditions, one might then ask if 
there is some connection between £(0) values for fermions and bosons with adjacent spins 
(taking local boundary conditions), reminiscent of supersymmetry. One would also like 
to understand whether the equality of the local and spectral values for £(0) is a feature 
peculiar to the highly symmetrical example of a 3-sphere surrounding a region of flat 
4-space, or whether there is an extension of this result to a more general context. 
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